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ABSTRACT 

We detected a major X-ray outburst from M82 with a duration of 79 days, an average flux of 
5 x KU 11 erg cm~ 2 s _1 in the 2-10 keV band, and strong variability. The X-ray spectrum remained hard 
throughout the outburst. We obtained a Chandra observation during the outburst that shows that 
the emission arises from the ultraluminous X-ray source X41.4+60. This source has an unabsorbed 
flux of (5.4 ± 0.2) x 10~ n erg cm -2 s _1 in the 0.3-8 keV band, equivalent to an isotropic luminosity of 
8.5 x 10 40 erg s _1 . The spectrum is adequately fitted with an absorbed power-law with a photon index 
of 1.55 ± 0.05. This photon index is very similar to the value of 1.61 ± 0.06 measured previously while 
the flux was (2.64 ± 0.14) x 10 -11 ergcm~ 2 s _1 . Thus, the source appears to remain in the hard state 
even at the highest flux levels observed. The X-ray spectral and timing data available for X41.4+60 
are consistent with the source being in a luminous hard state and a black hole mass in the range of 
one to a few thousand solar masses. 

Subject headings: black hole physics - galaxies: individual: M82 galaxies: stellar content - X-rays: 
galaxies - X-rays: black holes 



1. INTRODUCTION 

The bright X-ray sources in external galaxies, known 
as ultraluminous X-ray sources (ULXs), are of current 
interest because they represent either an unusual state 
of mass accretion, not or seldom seen from Galac- 
tic stellar-mass black hole X-ray binaries, or binary 
systems containing i nterm e diate-mass black hole s 
(Colbert fc Mushotzkyj |1999t iMakishima et all l2000f) . 
The ULX with the highest obser ved X-ray flux, CX OU 
J095550.2+694047 = X41.4+60 (iKaaret et al.ll200l . is 
located in the nearby starburst galaxy M82 at a distance 
of 3.63 Mpc. If the radiation from X41.4+60 is isotropic 
and from an accreting object, then the mass of the ac- 
cretor must be in excess of 500Af Q to avoid violating the 
Eddington limit. X41.4 +60 is offset from th e dynamical 
center of the galaxy (jKaaret et alJ l2001h . produces 
quasip eriodic oscillations (QPOs) in the 50-120 mHz 
range (IStrohmaver fc Mushotzkv 2003; M ucciarelli et al.l 
20061: iDewangan. Titarchuk. fc Griffiths! 120061 : 
Feng fc Kaaretl I2007T). produces little if any radio 



emission (|Kaaret. Simet. fc Lang] l2006bf), and exhibits 
an X-r ay periodicity of 62 days (jKaaret , Simet. fc Land 
l2006allbT Kaaret fc Fengl l2007h . The position of 
X41.4+ 60 is within 1" of the position of an infrared 
source (jKaaret et al.l 12004ft identified as the super star 
cluster MGG 11 (|McCradv. Gilbert, fc Grahan. l2003h 
that has an extremely dense core that may have led 
to numerous stellar collisions an d the formation of an 
interm ediate mass black hole (|Portegies Zwart et al.l 
l200l . 

In on-going monitoring of M82 using the Proportional 
Counter Array (PCA) on the Rossi X-Ray Timing Ex- 
plorer (RXTE), we detected a major flare from M82 with 
a peak flux corresponding to an isotropic luminosity of 
1.1 x 10 41 ergs~ 1 . After the onset of this remarkable X- 
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ray flare, we obtained imaging X-ray observations with 
the Chandra X-Ray Observatory. We describe X-ray ob- 
servations with RXTE in § 2, with Chandra in § 3, and 
discuss the results in § 4. 

2. RXTE OBSERVATIONS 

We monitored M82 usin g the PCA on RXTE 
(jBradt. Rothschild, fc Swankl I1993T ) under programs 
920 98 and 93123 (PI K aaret). As described 
in jKaaret. Simet. fc Lamp (|2006al lbl): Kaaret fc Fengl 
(|2007h . we analyzed the real time data from RXTE 
soon after the data became available on the RXTE ftp 
site which enabled us to request a target of opportu- 
nity observation with the Chandra X-ray Observatory 
(jWeisskopf et al.ll2002h after detection of the flare. 

The RXTE results presented here were obtained with 
the production data and the analysis was carried out 
with version 6.4.1 of the HEAsoft package and the 
CALDB available in June 2008. After the data were 
retrieved, they were filtered to select good time in- 
tervals such that Proportional Counter Unit (PCU) 2 
was on, the source was more than ten degrees above 
the horizon, the pointing offset from the source was 
less than 0.01°, the satellite was more than 30 s 
past the South Atlantic Anomaly (SAA), and the elec- 
tron contamination was low, specifically with an elec- 
tron rate less than 0.1 c/s. A background file was 
made for the same interval as the observation with the 
model pca.bkgd_cmfaintl7.eMv20051128.mdl. The tool 
saextrct was used to make both light curves and spec- 
tra for the source and background, using only data from 
PCU 2. Spectra in the 3-12 keV band were fitted with 
a power-law model using the spectral fitting program 
XSPEC to obtain a measurement of the flux in the 2- 
10 keV band. A power-law model with interstellar ab- 
sorption with the column density fixed to 1.2 x 10 22 cm -2 
was used. The average photon index for the observa- 
tions shown in Fig. [T] was 2.15 with a standard deviation 
of 0.21. The average uncertainty in the photon index 
was 0.20, so there is no evidence for spectral variation. 



2 



Kaaret, Feng, & Gorski 



Fluxes calculated with the photon index fixed to the av- 
erage value agree with those calculated from the fits with 
the photon index left free. We note that the PCA field 
of view is about 1° FWHM. Thus, the PCA spectra rep- 
resent the summed spectra of all the sources in M82 and 
not the spectrum of the X41.4+60 alone; the spectral 
index measured with the PCA is not expected to agree 
with that measured for X41.4+60 alone with Chandra. 
We have reduced the PCA fluxes by a factor of 1.18 be- 
cause the PCA calibration appears to give systematically 
higher flux es than previous and other current X-ray in- 
struments (jTomsick et al.lll999t ). 

Fig. Q] shows the observed flux in the 2-10 keV band 
versus time from spectral fits with the photon index fixed 
to 2.15. The X-ray flare begins at MJD 54247.2 when the 
flux rises above 5 x 10 -11 erg cm -2 s . The flux remains 
high, although highly variable, until MJD 54325.7. After 
that date, the flux returns to the usual range for M82 of 
~ 2 — 3 x 10" 11 erg cm" 2 s" 1 . The flare duration was 
79 days. The peak flux was 6.1 x 10 -11 erg cm" 2 s" 1 , 
the average flux was 4.6 x 10 -11 erg cm" 2 s _1 , and the 
root-mean-square (rms) variation in the flux was 9.7 x 
10" 12 erg cm" 2 s" 1 indicating an rms amplitude of 21%. 
The dura tion of the flare is s i milar to that of the flare 
found bv IRephaeli fe Gruberi (|2002ft . The level of vari- 
ability is higher, but that may, at least in part, be due 
to the limited coverage (only four observations) available 
to IRephaeli fc Gruberi (|2002D . 

Fig.[2]shows a hardness-intensity diagra m for M82. We 
chose energy bands similar to those used by lBelloni et al.l 
(2005) in order to facilitate comparison with previous 
results. Specifically, we used the count rate in the 3.6- 
21.6 keV band as the measure of intensity and the ra- 
tio of counts in the 6.5-11.0 keV band to those in the 
3.6-6.5 keV band as the measure of hardness. The plot 
includes those observations shown in Fig. [T] for which 
the hardness could be calculated to a fractional accu- 
racy of better than 0.15. The photons recorded by the 
PCA include diffuse emission and other X-ray sources in 
M82 that affect the intensity and hardness. However, the 
changes in intensity and hardness are likely dominated 
by X41.4+60. Fig. [2] indicates that X41.4+60 is usually 
in the hard state and remains in the hard state during the 
whole of the outburst. There is no significant change in 
hardness as the intensity changes over a factor of about 
4. 

We searched for rapid variability follow ing the proce- 
dures described in iKaaret fe Fengl (|2007h . No statisti- 
cally significant quasiperiodic oscillations (QPOs) were 
found. This may be because the new RXTE observa- 
tions are rather short, typically less than 1000 s, and 
have fewer Proportional Counter Units turned on com- 
pared to previous RXTE observations in which QPOs 
were detected. 

3. OBSERVATIONS AND RESULTS 

After det ection of a high X-ray fl ux level with the PCA, 
a Chandra (jWeisskopf et al.ll2002f l Target of Opportunity 
Observation (TOO) was triggered. The observation (Ob- 
sID 8190; PI Kaaret) began on 2007 June 2 at 02:41:29 
UT (MJD 54253.112). The observation was made using 
the Advanced CCD Imaging Spectrometer spectroscopy 
array (ACIS-S) and the High Resolution Mirror Assem- 
bly (HRMA). We obtained 52,768 s of useful exposure. 



The ACIS-S was used i n imagin g mode. As described in 
IKaaret. Simet. fe Land |2006b) for our previous TOO, 
we attempted to reduce pile-up by offsetting the target 
— 3.'57 along the Y-detector coordinate in order to image 
the source using a point spread function that spread the 
beam over many CCD pixels and operating only the S3 
chip in the 1 /8 sub-array mode to reduce the frame time 
to 0.441 s. We also re-analyzed the Chandra observation 
(ObsID 6097) of M82 that began on 2005 February 4 at 
23:34:40 UT (MJD 53405.982). 

The Chandra data were subjected to standard data 
processing and event screening (Ciao version 4.0.1 
using CALDB version 3.4.3). Neither observation 
showed any strong backgrou nd flares. Following 
IKaaret. Simet. fc Land (|2006b), 

we constructed images 
for the 2-8 keV band and used the wavdetect tool to lo- 
cate the X-ray sources. Fig. [3] shows the images and 
the source regions for X41.4+60 defined using the 2- 
8 keV band. Note that the sources appear elongated 
because they are off axis. The ellipses encompassing 
X41.4+60 have semi-major axes of 5.69" and 7.11", and 
semi-minor axes of 3.33" and 2.93", respectively for ob- 
servation 6097 and 8190. The background region is de- 
fined by a 8.49" by 3.81" box placed near the source, 
but avoiding other point sources and the readout streak 
for X41.4+60. In observation 8190, some photons from 
another bright source, X42.3+59, fall into the region 
defined for X41.4+60. When accumulating the spec- 
trum for X41.4+60, we excluded photons from the re- 
gion around X42.3+59, indicated by the dotted ellipse 
in Fig. [31 in order to minimize the contamination. Re- 
sponse matrices were calculated using mkacisrmf and the 
spectra were fitted using XSPEC 12 which is part of the 
HEAsoft package. 

We estimated the pileup in the spectrum for each ob- 
servation by summing the counts (with no energy se- 
lection) in a 3x3 pixel box around the pixel with the 
maximum number of counts and then calculating the 
counts per frame within the box. Observation 6097 
yielded 0.18 counts/frame and observation 8190 yielded 
0.31 counts/frame. Thus, within the central 3x3 pixel 
box, the fraction of piled-up events is 5.2% for observa- 
tion 6097 and 10.3% for observation 8190. The central 
3x3 pixel box contains 40% and 38% of the total events, 
so the overall pile-up fraction is about 2% and 4%, re- 
spectively for observation 6097 and 8190. 

Even though the overall pile-up is low, the central re- 
gions of the image in the second observation are strongly 
affected by pile-up. We chose to use two methods to 
model the pile-up. The first was to use the ACIS pile-up 
model in XSPEC version 12.4. Ox. We fitted the spectral 
data in the 0.3-8 keV band with an absorbed power law 
with and without the pileup model. The pileup model 
parameters were set to match the frame time of the Chan- 
dra observations. Since the observation is off-axis, the 
fraction of photons in the piled-up part of the spectrum 
is different from the 95% for an on-axis point source. We 
set this fraction (the parameter psffrac) to the ratio of 
counts, with no energy selection, in the 3x3 pixel box 
around the pixel with the maximum number of counts 
to the total counts in the source extraction region. The 
values are 0.40 and 0.38 for observations 6097 and 8190, 
respectively. The alpha values from the best fits are 0.33 
for observation 6097 and 0.25 for observation 8190. 
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Fig. 1. — Light curve of M82 obtained using the PCA on RXTE. The plot shows the flux in the 2-10 keV band calculated for each 
observation versus the observation date in MJD. The vertical dashed lines indicate the time of the Chandra observation. An X-ray outburst 
occurred from MJD 54247 to MJD 54326. 



TABLE 1 
X-Ray Spectral Fits 



(10 22 cm 



Flux 

(10 -11 ergcm~ 2 s _1 ) 



X 2 /DoF 



Observation 6097 - 


2005 February 4 








MARX input 


1.61 


1.15 


1.55 




MARX output 


1.65±0.03 


1.20±0.03 


1.27±0.05 


433.2/401 


" pileup 


1.73±0.05 


1.24±0.03 


1.30±0.06 


421.7/400 


Data 


1.63±0.03 


1.18±0.03 


1.29±0.05 


347.1/407 


" pileup 


1.67±0.05 


1.20±0.03 


1.30±0.07 


343.1/406 


Observation 8190 - 


2007 June 2 








MARX input 


1.55 


1.29 


3.20 




MARX output 


1.50±0.02 


1.32±0.02 


2.67±0.08 


596.6/456 


" pileup 


1.63±0.04 


1.39±0.03 


2.77±0.11 


533.9/455 


Data 


1.50±0.02 


1.31±0.02 


2.66±0.08 


587.8/463 


" pileup 


1.56±0.03 


1.34±0.03 


2.70±0.11 


575.2/462 



Note. — The tabic includes for both observations: the parameters input to the Marx, the best fitted parameters for fits of an absorbed power-law 
model and an absorbed power-law model with pile-up correction fitted to the output of the Marx simulation and to the data. The columns give: V 
- the photon index, Nh - the absorption column density, Flux - absorbed source flux in the 0.3—10 keV band; and the x 2 an d number of degrees 
of freedom in the fit. 



To check the modeling of the pileup, the observations 
were simulated using Marx version 4.3.0. The Marx pro- 
gram is designed to simulate a Chandra observation. We 
set up simulations matching the off-axis pointing and 
roll angle for each observation. We followed an itera- 
tive process to find input parameters to Marx that best 
matched the data: we assumed a set of input parame- 
ter for an absorbed power-law model (absorption column 
density, photon index, and normalization) , ran Marx us- 
ing those parameters, fitted the output from Marx us- 
ing an absorbed power-law model, and then compared 
those best fitted parameters to those obtained from the 
data. Because of the high source flux, the parameters 
extracted from fits to the Marx simulation without the 
pile-up model differ significantly from the input parame- 
ters due to the effects of pile-up even the off-axis pointing 
and the use of a l/8th sub-array. We then adjusted the 
input parameters to Marx and repeated the procedure 
until a good match was obtained between the fits to the 
output of the simulation and the data. 



The final set of parameters input to the Marx simula- 
tion, the results of the fits to output of that simulation, 
and the results of the fit to the data are presented in 
Table [TJ We take the input to the Marx simulations as 
our best estimate of the true source properties. Com- 
paring the Marx inputs with the fits to the data us- 
ing the power-law model with pile-up correction provide 
an indication of the level of systematic uncertainty in 
the modeling. Our best estimate of the photon index 
is 1.61 ± 0.06 for observation 6097 and 1.55 ± 0.05 for 
8190. The absorption-corrected flux in the 0.3-8 keV 
band in units of 10~ n erg cm~ 2 s _1 is 2.64±0.14 for 
the first observation, and 5.4±0.2 for the second. In 
the 2-10 keV band, the absorption-corrected fluxes are 
1.86±0.10 and 4.0±0.2 in the same units. Thus, the pho- 
ton index remained constant, within the uncertainties, 
while the flux chang ed by a factor of 2. We note that 
iBerghea et al.l ()2008[ ) reported an unusually hard spec- 
trum for X41.4+60 for an observation on 1999 December 
30, but this observation is strongly affected by pile- up 
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Fig. 2. — Hardness-intensity diagram of M82 obtained using the 
PCA on RXTE. The plot shows the count rate in the 3.6-21.6 keV 
band as a function of the hardness, defined as the ratio of counts in 
the 6.5-11.0 keV band to those in the 3.6-6.5 keV band. The plot 
includes those observations shown in Fig.[T]for which the hardness 
could be calculated to a fractional accuracy of better than 0.15. 
The source remained in the hard state for the whole outburst. 

with a pile fraction of at least 17% in the peak 3x3 pixel 
box. 

The spectra of some ULXs appear to h ave curvature at 
high energies, around 5 keV and above (iFeng fc Kaaretl 
I2005t IStobbart et all 120061: iRobertsI 120071 ), and it is of 
interest to search for spectral curvature from X41.4+60. 
Unfortunately, the PCA observations around the time of 
the Chandra observation are too short to make a useful 
measurement of the high energy extension of the spec- 
trum. Also, as noted above, the PCA field of view is 
larger than the angular size of M82, so the PCA spec- 
trum would not reflect that of X41.4+60 alone. 

We fitted the Chandra data for observation 8190 using 
various models with pile-up and interstellar absorption. 
We fixed the alpha parameter to the value found in the 
single power-law fit as this was necessary to prevent alpha 
from being driven to 1.0 (in which case all piled- up pho- 
tons are accepted as valid, negating the pile-up correc- 
tion). A Comptonization model, comptt in XSPEC, pro- 
vided no improvement in the fit over the power-law model 
and the best fit temperature was well above 8 keV. An 
exponentially cutoff power-law model, cutoff pi, gave a 
cutoff energy larger than 8 keV, i.e. above the Chandra 
band. A multicolor disk blackbody model provides a sub- 
stantially worse fit than the simple power-law model. A 
multicolor disk blackbody model in which the temper- 
ature varies with radius as r~ p , the diskpbb, provided 
an adequate fit with p = 0.57 ± 0.01 and an inner disk 
temperature of 6.5^0 g keV. However, difficulties with 
the physical interpretation of this model given the high 
apparent luminosity and high temper ature needed for 
X41.4 +60 have already been raised by iMivakawa et al.l 
( 2008) based on spectral fitting with Suzaku. They con- 
clude that the model is physically untenable and that 
the source appears to be in the p ower-law state. We find 
an even higher temperature than Mivaka wa et alj (|2008t ) 
which would exacerbate the problems. 

Fig.[5]shows light curves of X41.4+60 in the 0.3-2 keV 
and 2-8 keV bands for the 2007 observation. The same 
source region used to extract the spectrum was used to 
extract the light curve. The time bin size is 2048 times 
the CCD frame time, i.e. 903.2 s. No attempt was made 
to correct for pile-up. The source was variable on time 



scales of thousands of seconds in the 2-8 keV band, but 
remained relatively constant in the 0.3-2 keV band. We 
divided the data into high count rate and low count rate 
halves and fitted the spectra for each using an absorbed 
power-law model with pile-up correction. The absorp- 
tion column density is consistent, within errors, while the 
photon index for the high count rate part, T = 1.52±0.03, 
is somewhat harder than for the low count rate part, 
T = 1.63 ± 0.04. This would favor spectral pivoting as 
the origin of the variability. 

4. DISCUSSION 

The observed flux of X41.4+60 measured with Chan- 
dra on 2007 June 2 is 3.0xl0~ n erg cm" 2 s" 1 in the 2- 
10 keV band. The flux from the remainder of M82 is 
1.8x 10 -11 ergcm -2 s _1 . The flux of the whole galaxy 
as measured by the PCA just before the start of the 
Chandra observation was (5.0 ± 0.2) x 10 -11 erg cm -2 s -1 
which is in agreement. We adopt 1.8 xl0~ n erg cm -2 s -1 
as an estimate of the flux arising from sources other 
than X41.4+60 in M82 for the PCA observations. Then, 
the peak observed flux in the 2-10 keV band from 
X41.4+60 during the outburst is 4.3xl0~ n ergcm~ 2 s -1 
and the average flux is 2.8xl0 -11 erg cm -2 s -1 . Assum- 
ing isotropic emission from a source located in M82 and 
converting these absorbed fluxes into intrinsic fluxes us- 
ing the Chandra spectrum derived above, the peak lu- 
minosity in the 2-10 keV band is 7.6 x 10 40 ergs _1 , the 
average luminosity is 4.9 x 10 40 ergs -1 , and the total en- 
ergy released in the outburst is ~ 3 x 10 47 erg. Extending 
the energy band to the full band, 0.3-12 keV, detected 
from X41.4+60 during the outburst would increase these 
v alues by 70%. 

iRemillard fc McClintockl (|2006f ) have identified three 
main spectral/timing states of stellar-mass black hole X- 
ray binaries (BHXBs): the steep power-law state, the 
thermal dominant state, and the hard state. The Chan- 
dra spectra show that X41.4+60 is quite hard during 
both observations with a photon index near 1.6. While 
it is not possible to rule out spectral curvature above 
8 keV, the observed 0.3-8 keV Chandra spectra are in- 
consistent with those observed for BHXBs in the steep 
power-law state or thermal dominant states. The Chan- 
dra spectra of X41.4+60 are adequately described by a 
single power-law with a photon index of T = 1.6 and are 
consistent with identification of the source as being in the 
hard state - defined as having a spectrum dominated by 
a power-law component with 80% or more of the total 
flux and a photon index 1.4 < T < 2.1, and pronounced 
continuum timing noise with an integrated rms power 
greater than 10%. We note that the photon index at 
energies below a ny break or cutoff in the power-la w com- 
ponent is used bv lRemillard &: McClintockl (|2006f) for the 
state classifications. Thus, identification of X41.4+60 as 
in the hard state is robust, even if the spectrum shows 
curvature at higher energies. 

The spectrum remains hard when the flux increases 
to the brightest levels ever seen from the source. The 
Chandra data show that the source is variable dur- 
ing the second observation, but are inadequate to 
detect fast timing noise. An XMM-Newton obser- 
vation made on 2004 April 21 reveals broad band 
timing noise (iDewangan. Titarchuk, k, Griffiths! 120061 : 
iMucciarelli et al.l l2006f) which appears to arise from 
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Fig. 3. — X-ray images of M82 in the 2-8 keV band. The panel on the left shows the Chandra observation on 2005 February 4, while the 
panel on the right shows the new Chandra observation on 2007 June 2. The intensity scale is the same in the two images. In both panels, 
the solid ellipse shows the source extraction region, the rectangle is the background region, and the arrow points North and has a length of 
5". For the 2007 observation, the dashed ellipse shows the source region for the X-ray source X42.3+59. This region was excluded in the 
extraction of photons from X41.4+60. 




Energy (keV) 

Fig. 4. — Energy spectrum of X41.4+60. The top panel shows 
the counts cm~ 2 s _1 keV - 1 for both observations on the same 
scale. Observation 8190 is the upper curve at all energies. The 
middle panel shows the residuals of the fit for observation 6097. 
The bottom panel shows the residuals for observation 8190. The 
spectral shape at high energies is very similar between the two 
observations. The column density for observation 8190 is slightly 
higher than for observation 6097. 



X41.4+60 (|Feng fe Kaaretll2007t) . The integrated power 
in the 6-80 mHz range is near 17%. This value is con- 
sistent with the magnitude of the rms power required 
for a source to be in the hard state, but the frequency 
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Fig. 5. — Light curves of X41.4+60 in the 2-8 keV band (crosses) 
and the 0.3-2 keV band (diamonds) for the 2007 observation. The 
source was strongly variable in the 2—8 keV band. 

band is shifted relative to the 0.1-10 Hz band used by 
iRemillard fc McClintockl (f2006h . 

The PCA data show that X41.4+60 remains in the 
hard state during the whole of the outburst, see Fig [5J 
There is no significant change in hardness as the inten- 
sity changes over a factor of about 4. This lack of change 
in hardness is similar to the evolutio n of GX 339-4 in th e 
rising phase of its 2002/3 outburst (|Belloni et al.ll2005l ). 
This behavior suggests that the mass accretion rate in 
X41.4+60 never becomes high enough at any point dur- 
ing the outburst to push the system out of the hard state. 

Stellar-mass black hole X-ray binaries have been ob- 
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served to remain in the hard state at luminosities as 
high as L/L E ~ 0-3 for GX 339-4 (jZdziarski et all 
120041: iMivakawa e t all 120081) and as high as L/L E ~ 
0.1 for XTE J1550-56 4 (jRodriguez. Corbel, fc Tomsickl 
120031 : lYuan et al.ll2007l ) , where L is the source luminosity 
and Le is the Eddington luminosity. Above this lumi- 
nosity, the sources transition to a softer spectral state. 
Assuming that the same limit of L/L E < 0.3 applies to 
more massive black holes found in the hard state, then 
the 0.3-8 keV flux measured with Chandra would imply 
a mass in excess of 2000 M Q if the object is radiating 
isotropically. Converting the peak flux observed with 
RXTE to a 0.3-20 keV luminosity of 1.6 x 10 41 ergs" 1 
would imply a black hole mass in excess of 4000 M Q . Al- 
lowing for moderate beaming (note that a geometically 
thin accretion disk has a beaming factor of 2 when viewed 
on axis), could reduce these values by a factor of a few. 

Using the correlati on between X-r a y spe ctral shape 
and L/L E found by iShemmer et al.l (|2008f) for radio- 
quiet active galactic nuclei (AGN), the photon index of 
T = 1.6 would suggest that L/Le ~ 0.1. The highest 
luminosity measured with Chandra would then imply a 
mass of about 7000 Mq. However, mass estimation via 
r and Lx is only ac curate within a factor of about 3 
(jShemmer et al.ll2008[ ). 

Quasi-periodic oscillations have been de- 
tected from the ce nt ral re gion of M82 
(IStrohmaver fc Mushotzkvl 120031: IMucciarelli et all 
120061: iDewangan. Titarchuk. fc G riffiths 2006) and were 
localized to X41.4+60 (|Feng fc Kaaretl I2007D . The 
strength and coherence of the QPOs suggests they are 
of 'type C as defined for stellar-mass b lack hole X-ray 
binaries (|Remillard fc McClintock1l20060 . Type C QPOs 
occur in the hard and intermediate states. Thus, the 
detected QPOs are consistent with the assertion that 
X41.4+60 remains in the hard state. 

Estimation of a precise compact object mass using 
the QPO frequencies alone is uncertain since the QPO 
frequencies for individual stellar-mass black holes are 
known to vary by large factors, but the QPO fre- 
quencies can be used to place bounds on the allowed 
range of masses. Comparing the maximum QPO fre- 
quency of 0.113 Hz observed from X41.4+60 with XMM- 
Newton to the highest observed QPOs frequencies for 
GRO J1655-40, GRS 1915+10 5, and XTE J1550-564 
(|Remillard fc McClintock1l2006D by scaling mass linearly 
with QPO frequency, leads to an upper bound on the 
compact object mass of 2.7 x 10 4 M Q . If the QPO from 
X41.4+60 is identified as a low frequency QPO, then a 
similar comparison with the highest observed low fre- 
quen cy QPOs for the same stellar-m ass black hole bina- 
ries (jMcClintock fc Remillardl l2006f) suggests an upper 
bound of 1700M Q . We note that the ratio of the QPO 
frequency to the frequency of the break in the broad- 
band timing noise of X41.4+60 is consistent with that 
seen for low frequency QPOs. This, and the fact that 
strong high frequency QPOs are seen only in the steep 
power-law state of stellar mass black holes while low fre- 
quency QPOs are seen in the hard state, suggests that 
the QPOs from X41.4+60 are analogous to low frequency 
QPOs. 

Better mass constraints may be obtained by com- 
bining spectral and timing information. We note that 
several authors have used XMM-Newton data to infer 



the spectral shape of X41.4+60. This is unreliable 
since X42.3+59 contributes significant flux is not 
resolved from X41.4+60 with XMM-Newton. Thus, 
all previous work suggesting the presence of a thermal 
component in the spectrum or identification of the 
source as being in the steep power-law or intermediate 
state must be re- visited (IStrohmaver fc Mushotzkv 



20031 iFiorito fc Titarchuk! 12004 IMucciarelli et al 



2006t Ibewangan. Titarchuk fc Griffiths! 120061: 



Okaiima. Ebisawa. fc Kawaguchil 120061 : ICasella et al.l 
2008). Simultaneous timing and spectral measurements 
of X41.4+60 will require simultaneous observations with 
XMM-Newton and Chandra. 

Accreting black holes of all masses are known 
to produce radio emission while in the hard state. 
The radio emission, X-ray emission, and black 
hole mass are related via t he 'f undamental plane' 
([Merloni. Heinz, fc Pi Matted I2003D . As noted in 
iKaaret et al.l (|2001l ). there was a radio transient in 1981 
close to the position of X41.4+6 0, but no radio emission 
has bee n detected subsequently. Kaaret, Si met. fc Land 
(|2006bf ) place an upper limit on the radio flux of 0.7 mJy 
at 8.5 GHz on 2005 February 3, just before the Chan- 
dra observat ion giving a 2-10 keV lumi nosity of 3 x 
lO^ergs" 1 . iParagi. Garrett. Biggsl (|2006l) place an up- 
per limit of 67 /xJy at 1.6 GHz on 2005 October 27, but, 
unfortunately, there is no contemporaneous X-ray obser- 
vation. For upper limits on the radio fluxes, the 'funda- 
mental plane' relation can be re-cast to place an upper 
limit on the black hole mass: 



M < 2000M 







/ x 1.28 / r \ -0.77 / s 

\10Q fj,Jy J \^10 4 °org/sy \^3.6Mpcy 



2.56 



where M is the black hole mass, F is the flux at 5 GHz, 
Lx is the luminosity in the 2-10 keV band, and d is the 
distance. From the 2005 February data, the limit is M < 
10 4 M Q . Because there is no X-ray information available 
for 2005 October, any limit is less secure. Assuming an 
X-ray luminosity of 1 x 10 40 erg s -1 , would imply a limit 
of M < 12OOM0. The scatter in the fundamental plane 
is significant, about a factor of 10, so the limits from 
the above equation are uncertain by at least a factor 
of 3. However, simultaneous X-ray and sensitive radio 
observations may provide a means to constrain the mass 
ofX41.4+60. 

There are several other known ULXs with luminosi- 
ties near 10 41 ergs~ 1 (|Gao et al.ll2003l : IWolter et alJ2006t 
iMiniutti et al.ll2006t ) . All of these sources have hard spec- 
tra and are, together with X41.4+60, perhaps the best 
candid ates for X-ray detected intermediate-mass black 
holes (|Robertsl I2007D . Unfortunately, due to the large 
distances to most of these hyper luminous sources, de- 
tailed spectral information and information on any spec- 
tral changes correlated with luminosity is limited. This 
trend toward hard spectra appear s to continue down to 
luminosities around lO^ergs -1 (fFeng fc Kaaret] 120061 : 
ISoria et al.ll2007l : iBerghea et al.ll2008l ). For sources with 
multiple o bservations, particular ly X-16 and X-ll in the 
Ante nnae dFeng fc Kaaretl 120061) and the ULX in NGC 
1365 (ISoria et al.H2"ool 7 the sources appear to remain in 
the hard state even when the flux changes by factors of 
several. This suggests that these objects also remain in 
the hard state. This is evidence against the suggestion 



X-Ray Outburst from M82 



7 



that these sources are stellar-mass black hole X-ray bi- 
naries that enter an "ultraluminuous branch" of the very 
high or steep power-law state at high luminosities, since 
a transition out of this state would be expected at lower 
luminosities. 
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